Dynamical Coulomb blockade of multiple Andreev reflections 
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We analyze the dynamical Coulomb blockade of multiple Andreev reflections (MAR) in a super- 
conducting quantum point contact coupled to a macroscopic impedance. We find that at very low 
transmission the blockade scales as n 2 with n = Int(2A/eV), where V is the bias voltage and A is 
the superconducting gap, as it would correspond to the occurrence of shots of charge ne. For higher 
transmission the blockade is reduced both due to Pauli principle and to elastic renormalization of 
the MAR probability, and for certain voltage regions it may even become an antiblockade, i.e. the 
current is enhanced due to the coupling with the electromagnetic environment. 
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Introduction: Our understanding of coherent electron 
transport in superconducting nanostructures has experi- 
enced a remarkable progress in the last decade A 
unified picture of the dc and ac Josephson effects has 
emerged in which the concept of multiple Andreev reflec- 
tion plays a central role. In this respect, supercon- 
ducting atomic-size contacts have provided an ideal sys- 
tem in which theoretical predictions can be tested with 
a high degree of accuracy Q. A remarkable feature of 
the MAR mechanism in superconducting quantum point 
contacts (SQPC) is the occurrence of coherent transfer 
of multiple charge quanta ne in the subgap region with 
n increasing as 1/V |{J|. This is particularly striking for 
the case of atomic contacts where the current can then 
be carried by the coherent transfer of several electrons 
through a cross section of only few atoms. 

On the other hand, the importance of charging effects 
in nanoscale junctions has been stressed since the early 
90's 0. In these junctions with a very small capaci- 
tance quantum fluctuations result in a suppression of the 
current at low temperatures, a phenomenon known as 
dynamical Coulomb blockade (DCB) Q- This effect is 
strongly dependent on the effective impedance of the cir- 
cuit in which the nanoscale junction is embedded. One 
can expect these effects to scale with the charge of the 
carriers, which would result in a much stronger blockade 
when several electrons are transferred. A natural ques- 
tion that arises is to which extent this simple picture 
applies for the blockade of MAR in a SQPC. Moreover, 
one would like to determine the actual signatures of DCB 
that could be observed experimentally 

The aim of this Letter is to address these questions 
in the case of a single channel SQPC. For this purpose 
we have extended the microscopic theory of transport in 
the MAR regime to take into account charging effects 
and, more generally, the effect of quantum fluctuations 
in the phase introduced by the electromagnetic environ- 
ment. We show that while in the tunnel limit the on- 
set of MAR processes is strongly blocked (the relative 
blockade scales as n 2 ), as the transmission increases the 



blockade is progressively suppressed. Moreover, at high 
transmission one can even observe the opposite effect, i.e. 
an enhancement or antiblockade of the current in certain 
voltage regions due to the coupling with the environment. 

Dynamical Coulomb blockade in a normal QPC of ar- 
bitrary transmission was analyzed in Refs 0, E3 ■ I n 
Ref. |s| it was shown that one can establish a direct link 
between DCB and shot noise in this type of systems. The 
relation can be written in a compact way as [JjJ 
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dS 

dw'sign(w - u') — (u',V), 
(1) 

where < SI > denotes the correction in the current due 
to DCB, J(w) is the Fourier transform of the phase corre- 
lation function J(i) =< <f>(t)<f>(Q)) > — < (j> 2 >, related to 
the macroscopic impedance, Z(uj), characterizing the en- 
vironment by J(i) = 2e 2 /h J duReZ(uj){e luJt - l)/w and 
Si(uj, V) denotes the current noise spectrum of the QPC. 
In subsequent papers Kinderman et al. |ld| extended the 
analysis of interaction effects to the full current distribu- 
tion. In particular, in the second paper of Ref.[lfJ it was 
shown that at zero temperature and for a conductor with 
a linear current-voltage (IV) characteristic a relation of 
this type holds for the full distribution. However, as we 
discuss below, this relation does not hold exactly in the 
case of superconducting electrodes. The problem thus 
requires the use of more general methods. 

Theoretical formalism: As discussed in Ref. 0, a QPC 
coupled to an environment can be described by a Hamil- 
tonian of the form H = Hl + Hr + Ht + H env where 
Hl.r correspond to the uncoupled electrodes, which we 
describe as BCS superconductors, H env to the environ- 
ment taken as a collection of bosonic modes and 

Ht = X)cr v (c^CRo-Ae + h.c^j is a term coupling the 

leads. The translation operator A e = e 1 ^ , with <f> sat- 
isfying [Q, <j)\ — ie, takes into account the change in 
the charge Q of the environment associated with the 
transfer process. Within this model the normal trans- 
mission is given by r = 4(vW) 2 /(W 2 + v 2 ) 2 , where 
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1/Wir is the density of states on the normal leads at 
the Fermi energy while the current operator is given 

byi = f E.^LoUe-h.c.). 

The Keldysh formalism provides a general starting 
point to analyze charging effects not only on the mean 
current but on all the moments of the distribution. The 
main quantity to be determined is the cumulant gener- 
ating function S( X ) = ~ l n -?Xx) with 



where /;; ; •(/,./,. \: D? ,;/:;<;",. 
i>f_ i {t) being a matrix in Nambu space associated with 
the hopping from i to —i: 



Set,/} 



F( x ) =< ^exp j-~ J H T , x (t)dtVS c (0,t ) >, (2) 

where T c is the chronological ordering operator on the 
Keldysh contour and S c is the evolution operator of the 
system along the closed time loop. We take this contour 
to go from to to and back, to being the observation 
time (assumed to be much larger than any other char- 
acteristic time). In Eq.© the system is coupled to a 
"counting field" xMi which enters as an additional phase 
factor e lx (')/ 2 in Ht- This field changes sign on the two 
branches of the Keldysh contour, i.e. x(t±) — iX- The 
mean current and all the cumulants C n of the distribu- 
tion are obtained by repeated differentiation with respect 
to the X , i-e C n = -(-2ze)"pJ x=0 [13. 

As in Refs.0, ^} we assume that the external 
impedance is small (z = Z/(h/e 2 ) <C 1) and expand the 
generating function with respect to SHt — Ht — H^' , 
where H T °^ denotes the coupling term in the ideal voltage 
biased case. To first order in Z we obtain 



SS(x) = ~\ I dh [ dt 2 J(t u t 2 )K(t 1 ,t 2 ,x) 

^ J c J c 



(3) 



where K(t,t', X ) =< T c I x (t)i x (t')e * L H ^.x^ dt >. The 
correction to the mean current obtained by deriving SS 
with respect to x does not exactly coincide with Eq.lJIJ. 
In the normal state, where the energy dependence of the 
transmission can in general be neglected, both expres- 
sions are equivalent. However, in the superconducting 
case Eq.Q is only valid when Si(ui,V) is sufficiently 
smooth in the scale u>o of the typical energies characteriz- 
ing the environment, which in general requires « A. 

Applying Wick's theorem in Eq.ljSJ) we obtain 
an expression for SS in terms of the Keldysh- 
Nambu Green functions Gf'f (t,t' ,\) = —i < 



>, where ipj = [c'^^n 

and J a ^{t,t') =< f c (t>(t a )4>(t'p) > - < (j)' 2 >. The in- 
dexes i,j = ±1 denote the left and right electrodes (let 
L = 1 and R = —I). Then, Eq.@ can be written as 



ss( x )= ]r r 'dh r ' dhr-^hM) x 

a<f}=± Jo Jo 
J2 (i-i)Tr [Gtf(h,h,0)a z fff(h,h,x)a z ] (4) 
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(5) 



with v c l R (t) = (Vft L (t)) . By means of a double 
Fourier transform the Green functions can be expressed 
as G{<jj, to') — J2 n Go n (uj)S(oj — uj' + neV) [3j and thus 

SS( X )=toY J J ^ J duj'J a ^(uj-uj') x 

a, {3 

]T ( i -j)^Tr[^4.( W ,0)a z t^(^-x)^]- (6) 

±1 n 

To determine the Green functions entering in the eval- 
uation of SS(x) one has to solve the corresponding Dyson 
equation G = g + govoG, where the o product is a 
shorthand for summation over intermediate indexes and 
g are the unperturbed Keldysh Green functions of the 
uncoupled electrodes in equilibrium 



9ij = * 



ng a + (1 — n)g r 



n (g a — g r ) 



-ng' 



(1 - n)g a 



(7) 



where g a = g r '* = gl+fa x , with g 



(—uj/A)f, are the advanced and retarded BCS Green 
functions, and n(uj) is the Fermi factor. 

Results: In order to illustrate its effect on the IV char- 
acteristic we shall consider the simplest case of an envi- 
ronment characterized by a single mode of frequency loq. 
A more general situation can be straightforwardly ana- 
lyzed as a superposition of modes weighted by ReZ(w). 
Fig. 1 shows the overall behavior of the current blockade 
at zero temperature obtained by numerical evaluation of 
Eq.@ for ojo — 0.2 A and different values of the trans- 
mission. We also show for comparison the mean current 
< j > in the absence of environmental effects. The most 
prominent features are: i) in the tunnel limit the blockade 
appears around the threshold MAR voltages 2A/n, ii) as 
the transmission increases the region where blockade is 
effective moves towards V — ► 0, vanishing eventually for 
perfect transmission and iii) for high transmission and 
low bias the correction to the current < 51 > can ex- 
hibit a sign change. These features are discussed in more 
detail below. 

Let us first discuss the r — > regime. Here, at a fixed 
bias voltage V, SS is mainly determined by processes of 
order n ~ 2A/eV in the transmission. These processes 
are illustrated by the different diagrams in Fig. 2, which 
correspond to n = 3. Diagrams of type a) in Fig. 2 cor- 
respond to a real excitation of the environmental mode 
which gives a maximum contribution to SS in the volt- 
age range 2A/n < eV < (2A + u>o)/n, i.e. between the 
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FIG. 1: (color online) Mean current (upper panel) and cur- 
rent blockade (lower panel) for a single mode environment 
with loo = 0.2A at zero temperature for different transmission 
values. The current blockade curves have been displaced ver- 
tically for clarity. Notice normalization with Gn = (2e 2 //i)r 
and z, which is assumed to be a small parameter. 



opening of a MAR of order n until the onset of inelas- 
tic processes. More precisely, their contribution to 8S at 
zero temperature around eV ~ 2A/n is given by 



zen tv 



G) 



2d 



du>p(u>)p(u> + neV)K n (uj,uj) 



' A-neV 



-A— 
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with 



dujp(uj + u)q)p(u) + neV)K n (uj + loq, lo) > (8) 



K n (uj,uj') = 



e n + key ) n w + iev ^ 

j=l k=l l=n-j+l 



where p(oj) = lmg(uj)/TT is the BCS density of states. 
On the other hand, diagrams of type b) in Fig. 2 cor- 
respond to an elastic renormalization of the tunneling 
rates. In the tunnel limit they give a much smaller con- 
tribution which can be neglected for « A. From 
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FIG. 2: (color online) Typical diagrams contributing to the 
blockade in the tunnel limit eV ~ 2 A/3. Crosses indicate the 
hopping events on the two branches (+ or -) of the Keldysh 
contour. Full lines with an arrow indicate electron or hole 
propagators, while double arrowed lines correspond to the 
anomalous ones and wavy lines denote phase correlators. Di- 
agram a) is an inelastic correction due to the coupling with 
the environment while diagram b) corresponds to an elastic 
renormalization of the MAR probability. 



Eq. iJHJ and taking into account that in this limit < 

I>- TT (I) 2 n i A-neV duJ p(uj)p(co + neV)T n (uj), where 

r «M = Uk=i + keV )\ 2 0> we find that the Na- 
tive blockade — < SI > / z < I > reaches its maximum 
value n 2 for 2A/n < V < (2A + u )/n. This result is 
consistent with the fact that in the tunnel limit the trans- 
mitted charge is well defined and increases in a staircase 
way as the voltage is reduced 0. In this limit the en- 
vironment produces a blockade of MAR processes as if 
they would correspond to single " shots" of charge ne un- 
til their energy is sufficiently large to excite an environ- 
mental mode. 

There are, however, two basic mechanisms which re- 
duce the blockade at finite transmissions and finite ojq. 
First, elastic renormalization of the tunneling rates (di- 
agrams b) in Fig. 2) give a contribution to SS which in 
the tunnel limit and for small frequencies can be written 
as 



dujp(uj)p(uj + neV) 



I A-neV 



d 2 T n (cj) 
dcj 2 



(9) 



At the opening of the Andreev channels T„ exhibits a 
minimum as a function of energy and thus the contribu- 
tion to the current arising from SS^ is positive. On the 
other hand, as the transmission increases there is a sup- 
pression of blockade at large bias due to Pauli principle 
which is closely related to the reduction observed in shot 
noise 

The evolution of the relative blockade for r ranging 
from 0.01 to 0.4 is illustrated in Fig. |3| The maximum 
relative blockade expected in the tunnel limit is indicated 
by the shadowed regions. As can be observed deviations 
from the n 2 scaling are already significant for transmis- 
sions of the order of 0.1. 

The behavior of the blockade is qualitatively different 
for higher transmissions. This is illustrated in Fig. 
which shows < SI > close to the ballistic limit and for 
low bias voltage. It is observed that a blockade peak 
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FIG. 3: (color online) Relative blockade for r = 0.001, 0.1, 
0.2, and 0.4 (from top to bottom) and u>o = 0.2A. The shad- 
owed regions indicate the maximum relative blockade (scaling 
as n 2 ) which is predicted for r — » 0. 
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FIG. 4: (color online) Current blockade at small bias voltage 
and large transmission. Notice the sign change for interme- 
diate bias. The dotted lines correspond to the approximate 
expression of Eg. 1101 . The inset shows the behavior of block- 
ade for r = 0.9 and different values of Wo- 



at low bias is followed by a region where the system ex- 
hibits antiblockade. The overall behavior is robust with 
respect to changes in u>o from to ~ A (see inset in Fig. 
0J. One can understand this behavior from Eq.lJIJ which 
is valid when loq — > 0. In this limit and at zero tem- 
perature we have < SI >oc — uj 8 gy ^ ■ By virtue of this 
relation the blockade can be interpreted as a reaction of 
the system in order to compensate the power which is 
dissipated in the circuit by the current fluctuations in- 
troduced by the mesoscopic conductor. In a SQPC at 
high transmission the noise exhibits regions with nega- 
tive slope (dSi/dV < 0) 5]. In these regions the com- 
pensation requires an enhancement of the current instead 
of a reduction. This behavior can be obtained in a simple 
way as follows. Close to the ballistic limit the noise at 
low bias arises from Landau Zener transitions between 
Andreev states (see second Ref. in and is given by 
5(0) = ^Ap(l-p), where p = exp [-ttA(1 - r)/eV] is 
the Landau-Zener probability. As a result one obtains 



< SI >r 



2A 

^0 7/2 P 



(1-p) 



7rA(l 



eV 



(1 - 2p) 



(10) 

which exhibits a sign change for eV ~ 1.75A(1 — r). 

Conclusions: We have analyzed dynamical Coulomb 
blockade in a SQPC in series with a macroscopic 
impedance. A scaling with the square of the effective 
charge is found in the limit of vanishing transmission and 
environmental frequencies. As the transmission increases 
the blockade is strongly suppressed for eV ~ A and for 
certain voltage range it may even become an antiblock- 
ade. Our predictions could be tested experimentally in 
superconducting atomic contacts with a well character- 
ized electromagnetic environment, like those already used 
to analyze DCB in the normal state [l4j . 
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